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lo MATHEMATICS. CODING AMD APPLICATIONS 

1.1 Introduction 

During the period covered by this report 238 coded programs were mm. on 
the time . allocated to the Scientific and Engineering Computation (S&EC) Group^ 
These programs represent part of the work that has been carried on in 26 of the 
problems that have been accepted by the S&EC Group. Progress on 22 of these 
problems is given below in terras of programming hours, minutes of computer time, 
and progress reports as submitted by the programmers in question* 

No new problems were initiated during this period. Production runs are 
being carried out on seveh problems (#101, optical projperties of thin metal films? 
#109;, three-dimensional airplane pursuit-course; #112, solution of an eigenvector 
equation for the Educational Testing Service 1 #113* stress analysis of an L-t 
shaped homogeneous planar structure? #13^» determination of the eigenvalues and 
eigenvectors of symmetric matrices; #137, investigation of atmospheric turbulence 
by generalized harmonic analysis; #14?, study of energy bands in crystals). 

Professor F» Hildebrand of the MIT Mathematics Department spoke at the 
Seminar on Computing Machine Methods on November I?* His topic was "Optimum 
Interpolation*' - a group of about 75 interested persons attended* 

It is planned to begin the CS introductory programming course on December 
7* This is a resumption of a similar course that was given during the past 
year. A brief description of the course is given below under Problem #100. 

1.2 Programs and Computer Operation 

The following sxanmary is included as a guide for interpreting the abbreviations 
used below. A more detailed description of the terms involved can be found in 
M~2W* 

a« The upper case let^ter^ following the problem number has thi following 
significances 

A implies the problem is HOT for academic credit, is UNsponsored. 

B implies the problem IS for academic credit, is UNsponsored. 

C implies the problem is NOT for academic credit, ^ sponsored. 

D implies the problem IS for academic credit, IS sponsored. 
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The absence of a letter indicates that it is an internal S&EC problem* 

bo DIG denotes the Division of Industrial Cooperation* 
DCL denotes the Digital Computer Laboratroy* 
(MW denotes the Committee on Machine Methods of Canputatione 
DDL denotes the Division of Defense Laboratories* 



100 » Comprehensive System of Service Routines , developed by the S&EC Group 
at the Digital Computer Laboratory for the input conversion of 
suitably prepared punched paper tapes. When so requested, these 
routines automatically provide a program with suitable programmed 
arithmetic, cycle-counting, and output facilities* 
sDCL Staffs Ardenj 12 hours; Best, 40 hours 5 Demur jian, 19*5 hours; 
Denman, 10 hours j Frankovich, 20 hours? Hazel, 13 hours? Helwig, 70 
hours? Porter, 18 hours; Combelic, 15 hours; Wis 364 minutes 

A :few changes were made and tested to improve the form of the Floating 
Address (Flad) Table for Delayed Print programs. These changes will be incor- 
porated in the Comprehensive System during the next biweekly period* 

Two test programs were run to check four of the subroutines (viz*, print- 
out of WW instructions, interpretive instructions, decimal integers, octal 
integers) for the Drum Post«41ortem program* 

Hazel 

Further modifications for the Compreji:>ensive System have been tested* 
These changes which chfeck for some of the common, logical errors in the program 
being converted will be described in a forthcoming bulletin board memo* The 
title print program also is undergoing further changes to take into account the 
new form for tape titles* 

Frank ovich 

It is planned to begin the CS inti^oductory programming course on December 
7* This is a special two weeks course given by the S&EC Group and is intended 
to prepare selected personnel to program problems approved for solution on 
WWI, CS Manual I will be used as a text in this course <, The first eight 
chapters of this manual will be available for distribution by December 7* 

Denmanj Kopley, Porter 

Two test programs have been written to test the ability of the new PA 
to handle (24^6) numbers* The first uses every interpretive instruction except 
the new ics and prints out -results at each step as octal numbers* 

The second concentrates mainly on imr and idv, trying these orders for a 
large number of cases and printing out the cases with the largest round-^off 
errors* 

Best 

Work is continuing on the developnent of the new version of the compre- 
hensive system which will include a revised version of the (24,6) programmed 
arithmetic (PA) subroutine* It is planned to leave both the new version of CS 
and the old one on tape unit for a test period of at least one month* A 
brief bulletin board memo will be issued immediately defining the new PA* 

Helwig 
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101 C» Optical Properties of Thin Metal Films on transparent backings are 

determined and printed out automatically by this program ? the input 
data consist of the observed reflection and transmission coefficients s 
the index of the backings the wavelengthj and the sample thickness. 
The program calculates by means of an iterative procedure and prints 
out the rate of variation of these constants with reflection and 
transmission^ and the filja:i''s conductivity and dielectric constant* 
sfor Professor L, Harris,, Chemistry Department* Dr* A.L.Loeb 
sby Dr, A.L>Loeb,(DIC) , 15 hours g J» Richmond(DIC) jl5 hours 
sDCL Staffs Denman,, 2 hours » WIj 2? minutes 

The delayed printer modification and the column headings routine ncm 
operate successfully « A modificationp requiring only four registers, was pro- 
grammed and tested. This modification converts the R(ni,k)'='T(n5k) program into 
an n(R5T)»k(R5,T) program. This also worked satisfactorily using the delayed 
printers and it will be used to examine published data on fil:ms with thickness 
less than 75a for experimental accuracy* 

106 C. MIT Seismic Project is concerned with the development of methods for 

locating deep reflections frqri underground strata in seismic pros- 
pecting « The basic method is one of prediction by means of an 
optimum linear operator^, 
ttoT Professor P«M*Hurleyj Geology and Geophysics ^Professor Go Wadsworth, 
Mathematics Department 

°by Ej.A.»RobingQB(Res> Assoc*) » H. Briseoe> J6 hours; S» Smpson, 30 hours; 
Wo Walsh 8 20 hours 
iDGLi WWI, 255 minutes 

A new prediction program has been tested and used to form prediction 
errors for 50 linear operators. Thirty- two frequency spectra of seismograms 
have been computed in the past two weeks* 

Further tests of matrix inversion with a WWI subroutine were unsatls~ 
factory and a new programming approach has been started on this problem* 

10? C* (a) Autocorrelation and (b) Fo^xrier Transform^ Integral Evaluation * 
Programs were developed for these operations for the purpose of 
obtaining power spectra. The problem remains open for people who 
want to use these programs* 

sfor J.E.¥ardj Project Engineer pServomechanisTfis Laboratory 

•by DaX^Jsss(DIC), 30 hours; D* Hamilton, 5 hours 

sDCLs WWI9 k9 minutes 

During this period problem 10? was reopened for the purpose of locating 
and correcting the error in the Fourier Transform Program being used on problem 
136. At first it was thought that the errors were due to round^-off in the (158 0) 
recursion formula used to calculate cosine values* Therefore these calculations 
were changed to (3O9O) form^ a worthwhile change which is being retained* 

The error actually was due to the modification of the Standard Transform 
Program of problem 10? for use on problem I36* At' the time this modification 
was made it was assiamed that the subroutines of the program could be combined 
at will* Actually s, when only Simpoon's Rule is used^ the product of the function 
value times the cosine value must be made before the next cosine value is cal- 
culated, rather than after. This oversight has been corrected and the negative 
values of power no longer appear* The Standard Fourier Transform program was 
corrected in problem 10? o 



M=2540 page k 

108 C« An Interpretive Program is being deyeloped that will accept algebraic 

equations, differential equations ^ etc* expressed on Flexowriter 
Punched paper tape in ordinary mathematical notation (within 
certain limits imposed by the Flexowriter) as input and automatically 
provide the desired solution, 

sfor Dr., J,H„LaningoJroi, Instrumentation Laboratory 

sby J<.H.LaningoJr.(DIC)„N. Zierler(DIC).4Q hours 

;DCL Staffs Hazel, lo5 hours § WWI, 63 minutes 

About ?C^ of the time used this period was devoted to a final checkout 
of the main program^ The remainder of the time was spent in debugging the 
differential,»eqxiati on-solving subroutine* It is hoped that this subroutine 
will be made operative in the next biweekly period* 

109 Ct> An Airplane Pur suit°C curse Program is being developed which will take 

acco\mt of airplane dynamics and projectile ballistics and thus 
determine an airplane pursuit course in three dimensions « The 
problem consists essentially of solving 14 simultaneous non-linear 
differential equations by the Runge-Kutta Method which is of f ourth- 
order accuracy* 

2 for 14r« J,B*Feldraanj Instrwrientation Laboratory 

iby M«HeHellman(DIG) J, 10 hours 

sDCL"^"S^fTr''"Tfankovichj, 1 hour; Porter, 2 hours '^ Wil^ 95 minutes 

A successful run for the three-dimensional pursuit course problem was 
accomplished. A couple of preliminary test runs of about ? minutes each were 
made in order to check out the program before attempting -a complete run. 
The complete run required about 75 minutes and covered 5^ steps* each step 
consisted of a 10-foot increment in the Runge-Kutta Methods The results 
appear satisfactory particularly in regard to the closing of the slant plane 
triangle and the gravity drop triangle ttiroughout the run* 

It is planned to make a two-dimensional pursuit course run using initial 
conditions comparable to those used in the three-dimensional pursuit course run^ 
This will give a basis of comparison for the accuracy of the solution of the 
two-dimensional pursuit course method particularly in regard to time of flight* 

112 C « Lawlev' s -HQthod of Factor Analysis 4s sipplied to a correlation matrix 

obtained from psychological :t8st.3 and grades given at the Naval 
vAcademyo This requires the solution of a modified eigenvector 
equation for this matrix and is accomplished by an iteration of 
Hotelling" s method for solving an eigenvector pr obi erne 

sfor Dr^ ,F« M<, Lordj Educational Testing Service,, Princeton^ N<,J« 

sDCL Staffs Deimans 3 hours? WWIg 27 minutes 

The matrix problem was run again with assumed ranks of 9 and 10 for the 
correlation matrix* No alarms occurred and these results are now being examined 
at the Educational Testing Service* 

113 Go A Stress Analysis of an L'^^shaped Homogeneous Planar Structure is being 

made for the ease of a concentrated static load» Tnis structure 
is approximated by a framework of bars which will deform in the 
same manner as the prototype. This framework is then analyzed 
using the principles of virtual work and Southwell relaxation 
techniques,, Boundary conditions 'have been specified for the edge 
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of the framework so that the deformations of the model will cson<- 
form to the actual deformations of the strueture* 

sfor Professor J <,S .Archer „ Department of Civil and Sanitary Engineering 

«t)y So Sydney(Res<. Assist* GMMG), 60 hoiars 

sDCL Staffs Kopley^, 2 hours i WWIj, 435 minutes 

Production runs have been started on one section of the program. Inter- 
mediate results have) been obtained in the form of $^6 punched tapes which can 
be used directly in subsequent computations • Testing is proceeding on the 
remaining sections of the program. 

119 C« Spherical Wave Propagation produced by the sudden release of a spherical 

distribution of compressed air in the atmosphere is being studied 
by numerical means. This involves replacing a set of non-linear 
hyperbolic partial differential equations in 2 independent and 2 
dependent variables by a set of difference equations written along 
characteristics. An iterative procedure is used to solve these 
equations. 

sfor Professor C.G^Linj> Mathematics Department 

?^y A. Ralston (GMMC). 5 hours 

sDCL? ma^ 48 minutes 

The only run in the past two weeks produced a set of results which 
checked previous computation. 

Most computer work is being suspended for a couple of weeks while a 
paper on previous results is prepared^ 

120 D. The Aero thermopressor =■ This problem is concerned with the development 

of a device :^ or increasing the stagnation pressure of a hotg high 
velocity gas stream by means of evaporative cooling. The analytical 
investigation being carried %ut on liHiirlwind I involves the step- 
by-step solution of seven simultaneous non-linear differential 
equations which describe the thermodynamic and dynamic behavior ©f 
the compressible flow within the Aerothermopressor. 

sfor Professor A. H, Shapir©^ Department of Mechanical Engineering 

^^y B^D.Gavril (Ees. Assists) ;, X5 hours 

sDGLs WIj, 279 minutes 

Since last August ^ when computational work on this problem was tempor- 
arily suspendedj, a report on this investigation in the form of an Sc.D. thesis 
for the Department of Mechanical Engineering has been in preparation. This 
thesis 9 which is nearly completed 9 has indicated that significant improvement 
in aerothermopressor performance can be achieved by proper control of the Mach 
Number during the process. In order to substantiate these conclusions » the 
computational work was resiamed during the current biweekly period and^with 
the exception of one pending run,, has been successfully completed. A resume 
of this work follows. 

The original program^ written last April 5 was modified to include 
provision for the Flexowriter code in digits through 5 ©f the accumulators 
and several operating conveniences were added. These include a programmed 
stop character for the delayed printer and a carriage return counter for pro=^ 
viding proper spacing on consecutive sheets during typing. The various drum 
^^roll back*^ routines were retained,, and the 556 program tape^ which requires 
about 75 seconds for read in^ is automatically recorded on group 3 ©f the 
drum immediately following read in© 
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Two subroutines for controlling the variation of the Maeh lumber were 
written, tested^ and used for production runs. In this connections one of ^ 
these subroutines will be briefly described here to illustrate the flexibility 
and convenience of the programming procedure for Whirlwind I» 

In the original program, the simultaneous solutioh of the system of 
differential equations lead^ to a single governing equation of the form 

dbT "" 1«M^ ^^S'^dx^ ^-^^ 

where f and dg/dx are functions of the state of the stream at section x» The 
program was written only for computing dM^/dx from Eq, 1 for arbitrary function- 
al relationships A = A(x) and, in particular, for constant area ducts in which 
dA/dx = 0« In the calculations of the current biweekly, however, it was desired 
to study processes for which the Mach Number is varied in a prescribed manner, 
given 1:^ M = M(x)« Equation (1) then determines the area as a function of x, 
A ^ A(x)« To accomplish this calculation without a major and risky revision 
of the programs a subroutine based on the following rearrangement of Eq« (1) 
above was useds 

dM^ _ M^ f ^ + M^ dg 
dx 1«M^ dx ' l5P dx 



""HP ^ dx ^dx ^ 



2 2 
dM^ / dM 



dA/dx=0 



* dA „ dx ^ dx ^ 



» « 



The original program is followed in every respect, and the derivative 
(dM2/dx)j. /^j^Q is calculated from Equation(l) in the usual manner • At 

^ ' this point the subroutine is entered, the correct dM^/dx 
computed from the prescribed relation M=M(x) , and the area variation then' 
computed from Equation (2) above. The main program is then reentered, and the 
calculations proceed in the customary way. 

This subroutine, and a second one which prescribes the temperature 
difference between the liquid and gas phase as a function of x, were written 
as parameters and provide a convenient and valuable alteration in the cal« 
culation procedure. 

Various functional relationships of the form M=tl(x) were used* The 
results verified the prediction of improved performance and will be described 
in the final report. One additional case for a special type of function M^(x) 
is to be computed during the next biweekly j though no additional pr (^ramming 
will be done tin til the completion of the aforementioned report, 

123 C» Earth Resistivity measurements are used to calculate the Slichter 

kernel function which, in special cases, can be analyzed to give 
the actual distribution of resistivity* The method involves 
least-square fitting a set of polynomials to the measured surface='r 
potential function and integrating the product of this set and the 
zero-order Bessel function. 
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sfor P.M,i Hurley, Department of Geology and Geophysics, DIG 5<-6915 
'^ K. Vozoff (Res^ Assist,CMMC),30 hours 
sDCLs WIj 66 minutes 

The integration program, T 269^^ has been stopping on the calculation 
of the last of six polynomial coefficients which are used for interpolation 
on the normalized potential (r,p) curve. The reason is still undetermined. 

T 337? » the subprogram calculating JqCAt), is being run to check 
errors at various values of^ (X^)* Using a iSi of .1, it is in error lay .003 
at \ r*=ll if the calculation is started at ^ r=2. A 4 of .01 gives 5 place 
accuracy (at least) starting frcan X r^^O. 

126 G. A Data Reduction Program for use in the Servomechanj sms Laboratory is 
being developed in separate stage® to be combined at a later date. 
The first stage is concerned with devising a program to fit poly~ 
nomials to arbitrary empirical functions using a least squared 
error criterion. The procedure makes use of Legendre polynomials 
and matrix multiplication, 
sfor J. E.Ward, Servomechanisms Laboratory ,DIC No.7138,AF33( 616)2038 
rbv D.T.Ross(DIC).20 hours j D. Hamilton, 60 hours; R* Turyn,80 hours 
sDGL Staff s "Combelic, 2 hotirsi WWIj 123 minutes 

The Mistake Diagnosis Routine (MDR) was used to test sections of the 
Polynomial Fit Program, as previously described. These tests showed that the 
higher mofiients calculated by the Polytloiaial Fit Program diverge because of the 
fact that the buffer orders of the present Programmed Arithmetic (PA) Routines 
operate in (30»6) rather than (30,15)* This results in "large" zeros Which 
cause the divergence* The new PA being developed by Helwig will be substituted 
when it is ready^ 

The Lagrange Interpolation and Post~Mort@m programs have failed to 
work. Testing and modifications continue. 

The MDR report has been modified to refer to the new PA routines, but 
tapes are available for use with the present PA. Tape listings will be left 
in the Tape Room when the report is released in about two weeks. 

132 C. Subroutines for the Mumerical lv Contro] 1 fid Mini ngr Machine are being 

revised and tested. The set of subroutines facilitates programming 
I of the computations involved in the preparation of numerical data 

used to control the milling machine. The subroutines involve routine 
numerical and logical operations, 
sfor J. G.McDonough9 Sei^omechanisras Laboratory, DIG Project 6873 
sby J.H.Runvon (E«E.Res. Assist.), 20 hours 
sDGLs WI, 20 minutes 

The subroutine for finding series 16 camber line coordinates and 
slopes was successfully tested. Previous difficulties ,were caused by failure 
to specify a preset parameter in a library subroutine. ' 

Two more subroutines were written - a routine for selecting cut lengths 
for plane cuirves such that constant tolerance is obtained and a routine for 
evaluating the first and second derivatives on a series 16 wing section. The 
latter routine is used by the constant tolerance routine. A program for 
computing points on asymmetrical series 16 sections was also prepared and is 
being tested* 
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13^. C» Numerical Diagonalization Procedure > This program computes the exgen« 
values and eigenvectors of a symmetric matrix by a method of suc- 
cessive rotations. The program is available for use in any problem 
in which this calculation is required. 

2 for Professor J, C .Slater, Physics Department 

sby A. Meckler(DIC) 

sused by H. Kendall 

sDCL:Staffj Arden, 10.5 hours; Kopley,2 hours; Porter, ,5 hours 

Time was lost due to failure to write generalized decimal numbers 
correctly, a conversion error, an unexplained error in scope output, and the 
failure of the erase routine. 

After eliminating these difficulties, the eigenvalues and eigenvectors 
of five l^th order matrices were successfully calculated* 

137 D* Investigation of Atmospheric Turbulence as a noise input to airborne 

control systems. A stationary random process is assumed so that 
the methods of generalized harmonic analysis may be used to describe 
the turbulence components in terms of their power spectral densities. 

sfor Professor R.C.SearaAns, Department of Aeronautical Engineering 

sby R.A.Summers ( Res .Assist.). 5 hours 

sDCLi WWI, 55 minutes 

The Fourier transfonn program of problem 10? has been successfully 
corrected and several performance requests submitted. It is anticipated that 
all required transforms can be performed in about 2 hours of computer time. 
This should complete problem #137. 

138 B« Spheroidal Wave Functions are solutions of the scalar Helmholtz equation 

separated in spheroidal coordinates* A program has been developed 
for tabulating both the coefficients obtained by expanding the angular 
solutions of the first kind in associated Legendre functions and the 
coefficients obtained by expanding the 'radial solutions of the first 
kind in spherical Bessel functions. By analytically substituting 
these expansions in the separated ordinary differential equations, 
3- term recursion relations are obtained for the coefficients, and 
the radial and angular coefficients are found to be simply related. 
Both sets of coefficients are then determined by applying an iterative 
procedure to a continued fraction equation derived from one of the 
3-term recursion relations. The iteration proceeds until a value 
of the unknown separation constant of the differential equation 
which appears in the continued fraction equation is found such that 
it makes the coefficients compatible with their boundary conditions. 

sfor Professor P.MJlorse, Physics Department 

ibv F.J.Corbat^ CRes. Assist.CMMC),30 hours;J.D.C.Little(Res.Assist.CMMC) , 
30 hours 

sDCL Staff s Gombelic, 1 hour; IfWI, 8 minutes 

The final program has been partially revised in anticipation of the 
new GS PA being written by Helwig. An octal delayed print routine has been 
written for testing purposes* 
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IkQc Bmmev Session System consists of a conversion program 5 an interpretive 
routine, and mistake diagnostic routines stored in WI„ A special 
mnemonic instruction code has been developed for use with this 
system thus simulating a computer with characteristics quite different 
from those of WWI« This Summer Session (SS) computer was developed 
for the use of students participating in the MIT 1953 summer session 
course on "Digital Computers and Their Applications" ♦ The SS com- 
puter is being used in the E*E*Department courses 6»535 and 6,25 
and is available to programmers with suitable problems « 
sDCL Staffs Best, 37 hours? Combelic/ .19 hoiira Trankovichj 6 hours? 
Hoyo 38*5 hours? Siegel, 43 hours? WWI, 223 minutes 

The subroutines for adding a title-display feature to the SS computer 
have been written and are undergoing final tests* Only minor mistakes remain 
in the programs, and it is expected that these errors will be corrected and 
the final version of the routines incorporated into the fourth pass of the 
SS conversion program before the end of the next biweekly period* 

A temporary addition to the fourth pass, for the use of students in 
course 6«25 onlyj is being written* This subroutine will automatically produce 
identifying information and a set of calibrated co»ordinate axes on the 
oscilloscope superimposed on the curve to be plotted by the student's program* 



inated< 



A mistake in the instruction yin has been detected and will be elim= 



143 D« The Vibrational Frequency Spectihjm of a Copper Crystal is to be deter= 
• mined by solving a 3 x 3 secular determinant » each term of which 
consists of a finite Fourier Series -of 12 terms « This equation must be 
solved for 24,^5 different values of the wave-propagation vector* 
sfor Professors B.*E,,Warren and J,C*Slaterp Physics Department 
sby E*H^Jacobsen(Res* Assist J .1 hour? J«D,CoLittle(Res„Assist*CMMG),2 houri 
sD'CL I Wis 56 minutes 

The double integration concerned with second order temperature diffuse 
scattering of X-rays from a thenaally vibrating crystal lattice has been 
successfully completed* 

Also a frequency spectrum for the thermally vibrating lattice of 
^-iron has been successfully run making use of a simplified program previously 
developed by John D*C*34ittle. 

Meanwhile, work is progressing on a high speed program for determining 
the frequency spectrum of a thermally vibrating copper crystal lattice using 
exact tensor equations rather than the central force approximation employed 
in Little's program* 

1^ C* Self -Consistent Molecular Orbitals are the optimum choices of linear 
combinations of atomic orbitals determined through a process 
described as^ a self -consistent field approximation* The numeri©al 
procedure involves matrix- vector multiplications, vector additions, 
and matrix diagonalization. 

sfor Professor J*C ♦Slater, Physics Department 

*^y I^3r. Ao Heckler (DDL), 2 hours 

iBCL Staffs Arden, 7^3 hours? Porter, 2 hours 1 WI» 28 minutes 
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An obsciare error in the computation of the total energy was found and 
corrected. The convergence behavior of the routine is now being examined » 

14? C« Energy Bands in Crystals are being studied by finding solutions of the 
corresponding second order linear differential equation satisfying 
boundary conditions at the origin. The solutions are found approx^ 
iniately by using the Gauss^Jackson formula for forward integration* 
The solutions and their first derivatives are to be combined in a 
sum, the weighting factors being functions of an independent para= 
meter. 

sfor Professor J.CSlater, Physics Department j DIG No, 6853 

sby Dry D.J^Howarth CDIC) . 20 hours 

jDcL siafis Arden, 9 hours? Denman, 1 hour? WWI, 123 minutes 

The routine to calculate the Spherical Bessel Functions has been tested, 
made error free, and has been found to give 7-figure accuracy throughout the 
range required* 

The complete program has been assembled, and production work has started. 
In the interests of short periods of production time, it has been decided to 
divide the computation into two parts. First the integrations are performed 
for a set of energies- E, and the results stored on the magnetic drum § the 
computing timej including output, for one E is about ^ minutes, A routine 
has been written which will enable the results stored on the drum to be punched 
out on tape in a form suitable for input during the second part of the problem. 
This second part involves the computatiqn of a function of two variables E, 
Eq, The large number of Eq required can then be divided into ranges to enable 
shorter single production runs to be used. In this manner, it will be possible 
to select for a given range of E the range of Eq in which interest lies^ thus 
saving considerable computation time, 

1^ G, Digital Methods of Detecti wf^ Siffml frnm Mn-rs^ are being investigated, 
A sequence of binary numbers will simulate the Message wherein 
regions of high density of ones are signal regions and those with 
low density of ones are noise regions, Various methods of detect- 
ing the change from one region to another, as well as the length 
and midpoint of the signal regions are^ being studied, 

sfor J,V,Harringtonj, Lincoln Laboratory 

5 by G,P,Dinneen , Lincoln Laboratory 

sDGL? WCE, 3k minutes 

The two sequential observer detectors ahd the success run detector 
were tried for the weak signal case and for two signal intervals. For the 
short signal interval, the success run detector is more sensitive while the 
sequential observer with a bias of +4 is more, sensitive for the longer signal 
intervals, A modified sequential observer will be tried during the next period* 

152 D, Diffusion in an Qxide^Coated Cathode is a program to calculate the 

effects of combined thermal and electrolytic diffusion that occur 
in an oxide«-coated cathode when current is caused to flow through 
the cathode, 
sforW, B, Nottingham, Physics Department, DIG No, 6345 
sby H,B, Frost (Res, Assist, E> E, Department) , 30 hours 
sDCL Staffs Denman, 2 hours? Porter, 2 hours 1 WH.^ 3? minutes 
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A new program using an improved third- order approximation to the 
boundary conditions and an improved short routine for square rooting has been 
written and tested^ Unofrtunately, various errors have prevented a successful 
run thus far. Three coding errors, two errors in conversion, and one tape 
preparation error have been detected, and all known malfunctions have been 
explained. 

1*3 Operating Statistics 

Computer Time 

The following indicates the distribution of "WWI time allocated to the 
S&EC Group. 



Programs 


^1 hours. 


53 minutes 


Conversion 


14 hours, 


32 minutes 


Magnetic Drum Test 




18 minutes 


Magnetic Tape Test 




40 minutes 


Scope Calibration 




54 minutes 


Total Time Used 


58 hours. 


17 minutes 


Total Time Assigned 


63 hours, 


50 minutes 


Usable Time, Percentage 


91.3?^ 




Number of Programs, 


238 





1*^ Summary of Tape Room Bulletin Board Memoranda (I. Hazel) 

(These memoes are intended to inform programmers of changes in coding pro<» 
eedure, WI facilities", etc*) 

Planned Changes For CS 

Besides the changes in the new (30«-j,j) PA of CS, other changes listed 
below will be incorporated in the near future, 

1. The delayed printer will be used as output for the title and float- 
ing address (flad) table information instead of the scope camera. The letters 
CS will automatically precede all title information print-outs* The conven- 
tional form for titles may be used as well as "TAPE 100-34-5 JOSEPH GRJIPS" . 
Assigned and unassigned flads will be printed out in separate tables s the, 
assigned flads will be in the form "a6='483", and the unassigned flads in the 
fom "g7 at 48,509"* Example? 



CS TAPE HO 3004 M? JOSEPH GRMPS 
unassigned flads 



assigned flads 



zlO at 32, y9 at 33. x8 at 3^,^3»204 

al«32 a2«=40 a?=45 a9==58 
cl=93 e2=109 c5=82 
t4=204 tll=309 
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2« CS will convert the new instruction isc j ( j an integer) to isc 2j + 
initial address of counter block • 

3. The computer will stop and give an alarm when attempting to convert 
programs thati 

(1) overlap the PA or output routines that are being used, 

(2) make flad address assignment when the current address 
indicator is indefinite (see B-516-2). 

(3) contain the special word "(XJTAL"^ 

k» The cycle count block will be selected by the use of iti, iat, ict> 
ISLi or isc . as well as by the use of c in the program* 

Please forward any suggestions to J.M.Frankovich, Barta Building, Rm« 212, 
extension 631» 



1*3 Library of Subroutines 

The following is a list of subroutines available in. the S&EG Group 
Library of Subroutines. The routines are now arranged according to title 
rather than by tape number* Mod numbers are not needed since the taper oom will 
automatically use the latest mod* 

(Zk^S) Direct Print Generalized Decimal Number 

" Sinh X, Cosh x 

" . Sin Xj Gos X 

" Generalized Decimal Number Scope Layout 

" " " " Delayed Printout 
(30-jj) " " " 
Single Length Decimal Integer Scope Layout 
Fourth Order Runge Kutta, One Step, 'n» Diff, Eq* 
Auxiliary Buffer 
Largest Eigenvalue 

Single Length Decimal Integer Delayed Printer 
Simultaneous Equations (Craig's Method) 
Format and Delayed Print of Generalized Humbers 
Arc Sin x 

1*6 S&EC Library Accessions 

The f olloKring two books were added to the S&EG Library during the past 
two weeks* 

Har ydbook of Mathematical Tables and Formulas «Burington,R*S*« (B-I5) 

Mathematical Tables from Handbook of Chemistry and Physics «HQdgman,C*D*„(B^l6) 



Tape 


#3313 


M 


3328 


1! 


3329 


H 


3330 (*3395) 


n 


3332 


a 


3323 


n 


3333 


n 


2756 


H 


3322 


It 


2745 


11 


2767 


n 


3X68 


M 


321*^ 


« 


262^ 


M 


3217 


« 


3281 



* Tape #3395 is the new experimental version of the square root routine which 
is still under test. It is considerably faster than #3330» 
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2. COMPUTER ENGINEERING 

2.1 V/WI System Operation 

2»11 Core Memory (!♦ Holmes, A« J .Roberts) 

On Monday, November 16, the logic of the block orders bi and bo was 
revised^ One timing problem was encountered that was temporarily corrected 
with cabling changes. A permanent solution will be placed into effect in 
the near future. 

Also on November 16 the equipment of the computer room proper was 
modified to include the magnetic drums in the parity-check system. The drum 
equipment will be ready on Monday, November 23. 

Low margins within the pulse generator and the parity register have 
been corrected. The pulse generator required extensive modifications. The 
parity-register difficulty was encountered in digit and temporarily corrected 
by d-c coupling the cathode-follower outputs. Plans are being made to d-c 
couple the flip-flop and cathode followers of the entire register. 

2.12 Au3ciliary Drum (K. E» McVicar) 



Operation of the auxiliary-drum system has been quite good during the past, 
biweekly period. We have had no reports from operators of ar^ equipment 
failure during this time* 

2*13 Magnetic Tape (E* ?• Farnsworth) 

DifficTilties experienced by the operators with the delayed- output equip- 
ment because of errors in some programs are still being blamed on the equipment ♦ 
A check list and trouble-report form is being made up to insure more precise 
logging of suspected magnetic-tape-system troubles. This will serve the double 
purpose of aiding the system technicians ' in pin-pointing the trouble and giving 

the magnetic-tape group sufficient information to analyze or reproduce the 
difficulty. 

A Magnicorder reel-drive panel has been adopted for use with Raytheon 
spline type and triangular-hub reels as well as standard NAB reels to permit 
convenient filling of reels and transferring of tapes. Machine-shop work done 
here for the foregoing was so successful that duplicate parts are now being 
^turned out to convert unit 2 to triangular=l3ub reels at a fraction of the time 
and cost quoted by Raytheon* 

The noise existing in the servo loop of unit which caused small random 
movement of the reels and slack take-up mechanism has now been shifted to unit 
2 by interchanging servo-amplifier chassis. The noise arnplit|ide is now smaller 
and does not occur until unit 2 has been running for several hours. The trouble 
will be easier to track .down now that it is in the least critical tape unit. 

The remaining cross talk which was caused by the MTC unit-select switching 
transient getting into the delayed-printout system, has been cleared up by 
introducing a 5Q*-inicrosecond tiiae constant into the unit-select line from lOS. 
Considerable effort was expended to keep this voltage out of the printout system 
rather than change the MTC switching circuitry, but slowing down the switching 
proved to be the only solution. 
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The occasional reading diff icnlty experienced with unit during this 
period was traced by marginal checking to the read-record head. Microscopic 
investigation of this head revealed contamination by foreigVi particles? the 
reading difficulty cleared up after the head was thoroughly cleaned. In 
addition the problem of head wear, regrinding, polishing, etc, and the effects 
of this work on the reading-head magnetic air gap are being studied. Micro- 
scopic investigation has already produced some significant results. 

The problem of head reliability. is. also being undertaken by B. Paine, who 
is investigating the apparently inadequate electrical insulation in the Raytheon 
read/ record heads and will take up the problem with the manufacturer, 

2«14 Typewriter and Paper Tape (L,H,Korcott) 

During the past two weeks we have had several canplaints that one of our 
FL typewriters fails to carriage return properly while printing from magnetic 
tape. We have not been able to produce this failure in the shop, however, 
when operating the typewriter either manually or from paper tape, Farnsworth 
and I will again try this typewriter on magnetic tape in an attempt to track 
down the cause of this trouble, 

2«2 Terminal Eouipment (R, H, Gould) 

Two new types of film; for use with the Fairchild gcope camera were tested. 
They both have better contrast than the film now in use but are so much slower 
that the frame- counter number is not - recorded. The new camera control designed 
by N,L»Daggett will cure this trouble, and sharper, clearer pictures will be 
the result, 

A method of having fixed-scope gain position and intensity of different 
values available on a rotary switch is being devised. One se^iting would be 
for numbers formed of spots and the other for the smaller, dimmer numbers 
formed by the character generator. There is some risk of introducing more 
noise onto .the deflection lines, but care will be taken to avoid it, 

In-out system troubles in the last two weeks have consisted fyi part of a 
loose video cable, power«wiring error, junction«-box wiring error, and three 
timing troubles introduced by changes associated with the drum parity check « 
Equipment failures have been few* . 

2,21 Ferranti PETR (J,P,Stirman) 

The final |3readboard model of the Ferranti Photoelectric Tape Reader has 
been completed and tested. Its final test in ]ML will be delayed until the 
associated equipment can be constructed and the existing circuits modified. 

The tests consisted of plotting the output levels of each channel with 
the tube voltage as a variable. In this way the marginal performance of the 
circuit was checked. In addition, the output waveforms were observed on an 
oscilloscopej, with a tape input to the reader. The waveforms indicate that 
the per,formance of the unit will be satisfactory although there is no data 
available as yet on the effects of crosstalk between chafa^els* 
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^« ADMINISTRATION AND PERSONNEL 

Terminated Staff ( J. C, Proctor) 

Jerome Dint enf ass 
Peter Stephan ' 
Robert Garrett 

New Non~Staff (R. A. Osborne) 

James Ahlgren is an MIT student working part time in the Transistor 
Section of Group 62, 

Andrew Bowen has returned again to wrok in the Memory Section of Group 
62 on a part-time basis. 

Bessie Cachauni is a new messenger girl in the ^ittemore Building* 

Archie ,Lsmieux has joined the Inspection Department, 

Lucy MacFar land has returned to wprk in the Print Room, 

Thomas Murphy is a new member of the Construction Shop, 

Anthony Musi is a new member of the Inspection Department* 

Terminated Non-Staff (R .A •Osborne) 

Mildren Clark 
Gordon Morse 



